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In this article we discuss a model describing key features concerning the lineshapes and the 
coherent absorption conditions in Fano-resonant dissipative coupled oscillators. The model treats 
on the same footing the weak and strong coupling regimes, and includes the critical coupling concept, 
which is of great relevance in numerous applications; in addition, the role of asymmetry is thoroughly 
analyzed. Due to the wide generality of the model, which can be adapted to various frameworks like 
nanophotonics, plasmonics, and optomechanics, we envisage that the analytical formulas presented 
here will be crucial to effectively design devices and to interpret experimental results. 


Among the resonance lineshapes, Fano lineshapes de¬ 
serve a special interest since they first enabled to inter¬ 
pret atomic and molecular physical processes [1], and 
subsequently the response of nanostructured systems like 
photonic crystals and plasmonic resonators [2]. In gen¬ 
eral, resonances may be absorptive, in the sense that 
part of the energy vehicled by the excitation field is con¬ 
verted into other forms. While absorption, in the sense 
of losses, is usually an undesired effect, in other frame¬ 
works it can be harnessed to enable functional operations 
like detection, wavelength conversion or quantum state 
transfer. Absorbing systems also exhibit a rich physics, 
as for instance that of coherent perfect absorption (CPA) 
which shares some key mathematical aspects with non- 
Hermitian quantum systems and parity-time symmetry 
breaking [3-6]. In addition, CPA and related concepts 
may open new avenues in the control of wave properties 
like polarization [7] or in the processing of chaotic signals 
[ 8 ]. 

In its original formulation, the theory of CPA is very 
general, as it relies on fundamental properties of the scat¬ 
tering matrices [9]; when dealing with specific systems, 
appropriate modeling tools are needed. However, due to 
the possibly complex nature of the systems under analy¬ 
sis, microscopical approaches are often time-consuming, 
and an analytical model would be an advantage both for 
the interpretation of experimental results and as a guide 
to target ab-initio simulations. In this paper we analyze 
the CPA in a two-oscillator coupled-mode model, which, 
thanks to its generality, can be applied to a number 
of emergent frameworks like photon- or plasmon-exciton 
coupled systems [10-14] and optomechanics [15]. It will 
be first shown that the Fano trasmittance and reflectance 
resonances typical of a single symmetric resonator are in¬ 
herited by the coupled asymmetric system. Absorption 
lineshapes, instead, are described by another, universal 
lineshape, depending on few, physically meaningful pa¬ 
rameters. A check of the model validity is also provided. 


based on the concept of strong light-matter coupling in a 
realistic resonant metasurface embedding intersubband- 
active quantum wells. 



FIG. 1. Panel (a): sketch of the coupled oscillator model 
analyzed in the article. Panels (b-c): spectral lineshapes of 
the weakly (b) and strongly coupled (c) system. The trans¬ 
mittance lineshape can be tuned through a parameter (r, see 
text), and is inherited from the weakly to the strongly cou¬ 
pled case. The S-matrix determinant, instead, is given by 
a universal function, independent of the trasmittance (and 
reflectance) lineshapes. Parameter values are 7 r = 712 , 
LUc = u )\_2 = 507r, 7nr = 0. lu case (b) the ratio Tl/^r equals 
to 0.3, while in case (c) it equals to 2.7. 

The model under consideration is schematized in Fig. I 
(a). A resonant cavity at frequency Wc is coupled to a 
second resonant degree of freedom, here represented as a 
spring-mass resonator at through a coupling coeffi¬ 
cient H. From now on, the second oscillator will be re¬ 
ferred to as “matter” resonator, since a prototypical situ¬ 
ation would be that of a two-level system (atom, exciton) 
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treated under the semiclassical approximation. However, 
another important situation which can be described by 
the present model is that of a plasmonic system support¬ 
ing a “bright” and a “dark” (subradiant) mode: the os¬ 
cillator at uJi2 here represents the dark mode. The cavity 
resonator radiates into, and is excited from, two radiative 
scattering channels, with couplings Ki_ 2 - 712 describes an 
internal loss mechanism of the matter resonator, while 
7„r describes a non-radiative and non-resonant cavity 
loss mechanism. In the photonic framework, "fnr may 
represent losses such as roughness scattering or dissipa¬ 
tion in a metal component. 

The dynamics of the system is described by 

^ = (*wi2 - 712 )^ + i^a 

^ = {iujc - Jc)a i^b ((Kr)|s+) 

|s-) = (715+)-ka|«;) (1) 

where a and b are the amplitudes of, respectively, the cav¬ 
ity and the matter resonators. Here, |s+) = (5+,s+)^ 
describe the amplitudes of ingoing and outgoing scat¬ 
tered waves. Similarly, |k) = (ki, K 2 )'^■ In this notation, 
|n) —> (n| means transposition and complex conjugation, 
while |n) —)■ |n)* means only complex conjugation. C 
is a matrix which describes the non resonant scattering 
process. 

The free evolution of the cavity resonator occurs with 
a damping rate 7 c, which describes its total losses, ra¬ 
diative plus nonradiative. It can be decomposed as 
7c = 7r + 7 nr 5 where the second term has the mean¬ 
ing described above, while the first (the purely radiative 
damping rate) must satisfy 2jr = {k\k). This constraint 
can be derived by imposing instantaneous energy conser¬ 
vation to Eqs. 1. Energy conservation and time-reversal 
symmetry constraints also require C\k)* = — |k), as al¬ 
ready observed for the single-oscillator case [16]. 

The system’s response is fully described by its scat¬ 
tering matrix S', which links the amplitudes of ingoing 
and outgoing waves through jsg ) = S(w)|s[|"). Direct 
integration of Eqs. 1 yields 


S{uj) = C 


i{(jJ - U 12 ) + ^ 

(w — w+)(w — w_) 


( 2 ) 


where D = |/«) ((«:[*). The explicit expression of the poles 
u!± will be given in the following, while those of matrices 
C and D are given in the Supplementary Material. 

Matrix C depends on a parameter r € [0,1], which 
describes the off-resonant reflection amplitude. Matrix 
D involves an additional parameter which describes 
the asymmetry of the decay into the scattering channels. 
This parameter, constrained in the interval [—1,1], is con¬ 
nected with the coupling coefficients through the relation 
— l«^ 2 p)/(|KiP + |K 2 p)- This link between ki _2 
and r, already outlined for single-mode optical resonators 
[17], is here generalized to the two-oscillator model. 


In this model, if the cavity is decoupled from the mat¬ 
ter resonator (12 = 0 ) and there are zero non-resonant 
losses {"fnr = 0), the results of Ref. [17] are recovered. 
Considering transmission and reflection spectra, Fano 
lineshapes are obtained. Similar lineshapes, although less 
contrasted, are observed under the weak cavity-matter 
coupling regime, i.e., when the coupling D is finite but 
smaller than 712 and 7 c. As shown in Fig. 1 (b), the line- 
shapes sweep from a Lorentzian to an inverted Lorentzian 
according to the value of r. Here, we assumed ^ = 0, 
while it can be shown that ^ 7 ^ 0 leads to a further de- 
crased lineshape contrast. Consider now the spectra in 
Fig. 1 (c), obtained when O is larger than 712 and 7 c. 
The two resonators are strongly coupled, and the spec¬ 
tral feature is doubled, following the peak, dip, or asym¬ 
metric shape of the weakly coupled system. This unique 
behaviour {“lineshape inheritance"), which was already 
observed experimentally and justified heuristically [18], 
is now grounded on a basic theoretical model, and can 
be extended to all systems which can be described by 
Eqs. (1). 

In the above, the focus was on transmittance and re¬ 
flectance. While in an experiment these are the most 
easily accessible quantities, a more meaningful probe of 
a driven linear system would rather be the absorption, 
since it is directly connected to the excitation of the res¬ 
onant degrees of freedom. Here, the model under consid¬ 
eration presents two scattering channels, and should be 
analyzed in view of the coherent absorption theory for 
asymmetric two-port systems [9, 19]. A key quantity is 
the A-matrix determinant, which, following Eq. 2, can 
be expressed as det S{oj) = e^*‘^(a; — w+)(a; — w_)/(w — 
w+)(w — u!-). (j) is a global phase (see Supplementary 
Materials), while the zeroes, which are connected to the 
coherent perfect absorption (CPA), explicitly read 

Wc-fa;i 2 .7r--7«r-7l2 

^-*-X- 

2 2 

± ^ \J [(Wc - a;i 2 ) - f(7r - 7nr + 7l2)]^ + 

(3) 

The poles uj± are obtained from the zeroes by replacing 
7 r with — 7 r. The essential feature is that det S depends 
neither on r nor on the lineshape-governing factors 
do not influence the S'-matrix determinant. Nor they 
influence the CPA condition, too: it can be shown that 
det S' = 0 if and only if 

72 - 52 _ ^ y ^(^2 _ ^2 _ 4 f^ 2)2 + 4 ^ 2^2 

(4) 

where 7 . = 7 ^ - 7 „,. - 7 ^ 2 , 7 + = 7 ^ - 7«r + 7 i 2 , and 

6 = UJ12 — U!c- 

Eq. 4 unifies and generalizes the weak and strong crit¬ 
ical coupling concepts (WCC and SCC), which were in¬ 
troduced in Ref. [20] for a symmetric (^ = 0) and degen¬ 
erate ((5 = 0) coupled resonator system. As recalled in 
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Fig. 2(a), the WCC and SCC regimes appeared there to 
be well separated curves on the system’s phase diagram, 
merging into a single exceptional point. 



FIG. 2. Coherent perfect absorption (CPA) in the non¬ 
degenerate coupled oscillator model. The panels on the left 
represent the phase diagram of the system: on the thick 
curves, CPA occurs. Well separated weak and strong criti¬ 
cal coupling (WCC and SCC), encountered in the degenerate 
case (<5 = 0), are no more distinct in the S ^ 0 case. Panels 
(b) and (c) represent the S-matrix spectra for the parametric 
paths labeled B and C in panel (a), for a fixed value of S. 
Similarly, panels (e) and (f) correspond to points labeled E 
and F in panel (d), where instead 5 is swept. 

When instead i5 7 ^ 0, the solution of Eq. 4 are rep¬ 
resented by the thick lines in Fig. 2(a). Hence, for the 
non-degenerate case, well separated WCC and SCC do 
not exist anymore, and a smooth transition is instead 
observed. However, some features reminiscent of WCC 
and SCC are still evident when the det S spectra are 
considered. These spectra are reported in panels (b-c), 
which correspond to the parametric paths labeled B and 
C in panel (a). In case B, dots' has still a double-dip, 
but, as opposed to the <5 = 0 case, CPA never occurs 
simultaneously for two frequencies. When the transition 
between the SCC and WCC region is explored (path C), 
the coalescence of | det S | minima is observed, but only 
the higher frequency dip is a CPA. (Whether the CPA 
originates from the dip at higher frequency, or from that 
at lower, is determined by the signs of 6 and of 7 _). The 
effect of a continuous sweep of the detuning S is analyzed 
in Fig. 2 (d-f). Suppose that the system is close to, but 
not exactly on, the SCC (point E). Its spectrum, if J = 0, 
has two dips which do not reach zero, as highlighted by 
the dashed horizontal line in (e); if S is tuned, a CPA 
occurs either on the lower or on the higher frequency res¬ 
onance of the coupled system. Suppose instead that the 
system is close to the WCC (point F). Again, by tuning 


6 the system can be brought to CPA, but now with a 
single isolated spectral feature. In essence, the detuning 
has a twofold role. On one hand, it weakens the distinc¬ 
tion between strong and weak critical coupling, as when 
6^0 the curves describing CPA on the phase diagram 
are smooth and do not exhibit any exceptional point. On 
the other hand, a proper tuning of 6 can help a system to 
reach CPA, without acting either on the coupling or on 
the damping rates. All these observations apply indepen¬ 
dently of the system’s asymmetry and of the specific Fano 
lineshape observed in reflection/transmission, leading to 
a universal behaviour of coupled dissipative resonators. 

Nonetheless, the asymmetry plays a singular role in 
the response of coupled dissipative oscillators. This will 
be explored by analyzing a realistic structure, a res¬ 
onant metasurface embedding semiconductor quantum 
wells (QWs). Such device implements a prototypical sys¬ 
tem in which weak and strong coupling have been ob¬ 
served, and where they can be harnessed to develop effi¬ 
cient mid-infrared and terahertz light sources, as well as 
functional optical components [ 21 , 22 ]. 

A schematic of the device is sketched in Fig. 3 
(a). It consists of a heterostructured membrane of 
60 GaAs/Alo.aaGao.eyAs QWs, with well/barrier thick¬ 
nesses 6 . 8/20 nm resulting in an intersubband transition 
frequency ?iwi 2 — 150 meV. The membrane has thickness 
^2 = 1.3 fim, and is periodically patterned with 50 nm 
thick gold stripes, whose spacing is a ~ 3.5 /rm and filling 
fraction is / ~ 0.8. A high-index coating (e = 10) with 
thickness ti completes the stack. The structure is mod- 
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FIG. 3. Panel (a), schematic of a resonant metasurface 
embedding quantum wells, which implement a prototype of 
strongly and critically coupled oscillators with asymmetry. 
Panels (b) and (c), resonant field and spectral transmittance 
of two structures which differ by the value of ti. These res¬ 
onators are described by different asymmetry parameters, re¬ 
sulting in differently contrasted lineshapes. The calculations 
from a rigorous electromagnetic solver (dots) are faithfully 
reproduced by the coupled-mode model (red and blue lines). 


eled through rigorous coupled wave analysis (RGWA), 
following the details given in [18, 23] and in the Supple¬ 
mentary Material. If the response of the QWs is turned 
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off (ff oc V An = 0), the device exhibits an isolated 
photonic resonance at Wc — 150 meV (dotted traces in 
Fig. (3)). Two structures have been analyzed, and la¬ 
beled 1 and 2; their details are given in the Supplemen¬ 
tary Material. Their main difference is in the value of ti, 
which governs the cavity asymmetry. The RCWA trans¬ 
mittance spectra have been fitted through Eq. 2, which 
led to ~ 0 for struct. 1 and ~ 0.7 for struct. 2. A 
signature of such non-zero value is the transmittance con¬ 
trast, which, for struct. 2, is less than unity. When the 
QW response is turned on (An ^ 0), the transmittance 
peak splits into two, with the line shape and contrast be¬ 
ing inherited, ft should be noticed that, for the 0 
case, the agreement between the RCWA spectra (dots) 
and those obtained from the coupled mode model (lines) 
relies on a first-principle relation which connects An and 
n (see Supplementary Material for details). 

The most interesting physics occurs however when the 
effects of the asymmetry on the joint absorption Aj.± 
of the two-sided cavity are analyzed [19]. Absorption of 
optical radiation from such systems is a meaningful figure 
of merit, for instance for detectors and for other devices 
where an efficient pumping of the polariton population is 
needed. From Eq. 2, one gets 

Ay,± (w, x) = + r^x ± ^{1 - - x"^)^ . 

(5) 

Here, x = (|s]^P — P)/(|s]'’P -I- |sjp) is the asymme¬ 

try between the input intensities, and the sign ± rep¬ 
resents the minimum and maximum absorption achiev¬ 
able by acting on the relative phase of the input beams. 
Notice that Aj as a function of x shows the peculiar el¬ 
liptical behaviour of coherent absorption [19]. The sig¬ 
nificant feature standing out from Eq. 5 is that the el¬ 
lipse is described by the sole function in parentheses, 
which factors out from the universal spectral lineshape 
Auni(w) = 1 — |det«S'(w)p. Hence, it is only the func¬ 
tion in parentheses, which we label f±{r,^,x), that sum¬ 
marizes the effect of asymmetry on coherent absorption. 
Instead, Auni is independent on both r and This obser¬ 
vation goes beyond what has been stated when discussing 
Eq. 4: it is not only the CPA condition which does not 
involve r and but the joint absorption lineshape in its 
fullness. 

In Fig. 4 (a) and (b) we plot f±{x), where the choices 
= 0 and = 0.7 recall the actual values of struc¬ 
tures 1 and 2. The function /_ is zero for x = —r^, 
which means that the system can always exhibit coher¬ 
ent perfect transparency (CPT) provided that the excita¬ 
tion intensities properly match the intrinsic asymmetry 
of the photonic resonance. For what concerns CPA, in¬ 
stead, necessary and sufficient condition is that Auni = 1 
(and hence | dots'! = 0 )j f+ — 1) ^ 

other words, there are two independent requirements: the 
first on damping rates, detuning, and coupling coefficient 
(Eq. 4), the second on the symmetry. 





FIG. 4. Coherent absorption universal factor f± for a coupled- 
resonator system which behaves symmetrically (a) or asym¬ 
metrically (b). Minimum and maximum joint absorption 
for structure 1 (c) and structure 2 (d, e). Both structures 
show coherent perfect absorption and transparency (CPA and 
CPT), but in the latter those are observed at different x- 
values, i.e., for different states of external excitation, due to 
the asymmetric behaviour of the cavity. 


A numerical test is proposed in Fig. 4 (c-e), where Aj 
calculated numerically for structures 1 and 2 is compared 
with the prediction of Eq. 5. No further parameters in 
addition to those which led to the fitting of Fig. 3 are 
involved. Since the samples have been designed in order 
to have J = 0 and jr = 7 i 2 , in both cases strong critical 
coupling occurs and CPA is expected. Indeed, when the 
proper x-value is chosen, Aj reaches 1 in a double-peaked 
fashion for both structures (Fig. 4 (c) and (e)). Since 
structure 1 behaves symmetrically (^ = 0), CPA and 
CPT occur simultaneously for a; = 0; in structure 2, in¬ 
stead, CPA and CPT occur for opposite values ^ = ±0.7, 
consistently with Eq. 5. 

In summary, we studied the absorption lineshapes oc¬ 
curring at the transition between weak and strong critical 
coupling for a system consisting of two coupled detuned 
resonators, one of which is radiatively coupled with the 
exterior in an asymmetric fashion. From this model a 
peculiar fingerprint in the absorption spectra stands out: 
a universal lineshape, independent of the asymmetry de¬ 
gree, which instead rules the Fano lineshapes observed in 
transmittance or in reflectance. Similarly, the coherent 
perfect absorption (CPA) condition results to be inde¬ 
pendent of the lineshape-governing factors. Rather, it 
turns out that the CPA phase diagram is significantly 
affected by the oscillator’s detuning, whose effect is to 
weaken the distinction between weak and strong critical 
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coupling regimes. Being the present model of wide gen¬ 
erality, it is of significance in the development of many 
active research areas like nanophotonics, plasmonics or 
optomechanics, where the modeling of a complex sys¬ 
tem would draw advantage from analytical expressions 
depending on few parameters of direct interpretation. 
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These Supplementary Materials contains explicit expressions for certain quantities introduced in 
the main text, as well as detailed information about the analyzed nanostructured devices. 


EXPLICIT EXPRESSIONS OF MATRICES C 
AND D 


The coupled-mode model presented in the first part of 
the main article involves two matrices, C and D, which 
represent the non-resonant and the resonant scattering 
processes, respectively. The elements in these matri¬ 
ces are constrained by energy conservation, time-reversal 
symmetry and reciprocity. Under these hypotheses, the 
most general expressions explicitly read 


C = 


rC'l’ it \ 
it re-^^ ) 



where 

dll = (r + e ± it^l - 

d22= {r-^±it./l-e) 

di 2 = ±r\/l - + it (1) 

The phases (p and ip can be chosen at will, while r is 
a real parameter belonging to [0,1]. The parameter t is 
connected to r via t = -y/l — {t G [0,1]). The meaning 
and the bounds on ^ are discussed in the main article. 


METASURFACE MODELING AND 
PARAMETRIZATION 

The metasurface is modeled by means of rigorous cou¬ 
pled wave analysis. TM polarization (i.e., incident elec¬ 
tric field polarized orthogonal to the metal stripes) and 
normal incidence is considered. Permittivity of gold 
is assumed to be eau = —4000 -I- 300i, that of GaAs 
£GaAs = 10, while the quantum well’s (QWs) response 
follows a Lorentz oscillator model whose details are given 
in Refs. [1, 2]. There, the key parameters are the reso¬ 
nance frequency uji 2 , the damping rate 712 and the sub¬ 
band surface charge difference An. 


First we considered the QWs to be inactive (An = 0), 
which means that the metasurface resonance is not cou¬ 
pled to the QWs {ft = 0). We designed two structures, 
whose geometrical parameters a, / and ti are reported 
in Table I, as well as the parameters ojc, 7 r, r and ^ 
which describe the response in terms of coupled mode 
model. The thickness ti is the most significant struc¬ 
tural parameter, since it rules the transmittance contrast 
which is directly connected to the asymmetry parameter 
It turns out that ti = 0 leads to a fully contrasted 
transmittance lineshape, which corresponds to ^ = 0 and 
hence to equal decay rates in the two scattering channels 
(ki = K 2 in the notation of Fig. 1 of the main text). The 
parameters a and / have been finely tuned in order to 
obtain resonance frequencies and linewidths as close as 
possible for the two structures. By fitting the transmit¬ 
tance spectra of the uncoupled metasurface (empty dots 
in Fig. 3 in the main text), the parameters oJc, 7 r, e and 
^ reported in Table I are obtained. It can be noticed 
that both structures have similar radiative lifetimes 7 ^ 
and non-resonant reflection r. Moreover, both structures 
share similar values for the field overlap factors F. This 
is defined as F = / |Epe (QW stands here 

for the quantum well region), and will be of relevance in 
the following. 

In a second step we considered the active QWs (An = 
5 X 10^^ cm“^). In this case, the metasurface resonance 
is coupled with the QW transition, with a strength quan¬ 
tified by the vacuum Rabi frequency which reads [1] 


fi = 


e^AnF 


4EoEGaAs 


where £GaAs is the well material permittivity, m* = 
0.067mo is the conduction subband effective mass, and 
Lper is the QW period thickness. For the analyzed struc¬ 
tures, hit ~ 8 meV, and the system is in strong cou¬ 
pling. The fingerprint of strong coupling is evident from 
the double-peak transmittance curves obtained from the 
RCWA (dots in Fig.3, main text). In the RCWA the 
QWs are modeled as a Lorentz oscillator with huji 2 = 
150 meV, h'yi 2 = 3 meV. These numbers, together with 
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the other parameters extracted from the empty cavity 
spectra, are then plugged into the coupled-mode model, 
resulting in the transmittance curves reported as red lines 
in Fig. 3 of the main text. Such double-peaked spec¬ 
tra should be compared with the dotted traces obtained 
through the RCWA. The very good agreement witnesses 
that the analytical expressions provided by the coupled¬ 
mode model can quantitatively predict the response of 
an asymmetric two-port photonic system strongly cou¬ 
pled with a matter excitation. 


TABLE I. Geometric data and coupled-mode model parame¬ 
ters for the structures analyzed in Figs. 3 and 4. 



ti a f 

UJc r ^ 

r 

Struct. 

[pm] [pm] 

[meV] [meV] 


1 

0 3.70 0.80 

149.60 3.25 0.996 0.001 

0.47 

2 

0.5 3.20 0.73 

149.70 2.77 0.981 0.72 

0.56 
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